Extremal dynamics model on evolving networks 
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We investigate an extremal dynamics model of evolution with variable number of units. Due to 
addition and removal of the units, the topology of the network evolves and the network splits into 
several clusters. The activity is mostly concentrated in the largest cluster. The time dependence of 
the number of units exhibits intermittent structure. The self-organized criticality is manifested by a 
power-law distribution of forward avalanches, but two regimes with distinct exponents r = 1.98±0.04 
and r' = 1.65 ±0.05 are found. The distribution of extinction sizes obeys a power law with exponent 
2.32 ±0.05. 

PACS numbers: 05.40.-a, 87.10. ±e, 87.23.Kg 
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Extremal dynamics (ED) models |^ are used in wide 
area of problems, ranging from growth in disordered 
medium ||], dislocation movement friction to bi- 
ological evolution jsj. Among them, the Bak-Sneppen 
(BS) model [|j plays the role of a testing ground for var- 
ious analytical as well as numerical approaches (see for 
example JT|,|6[-|TTt ) . The dynamical system in question is 
composed of a large number of simple units, connected in 
a network. Each site of the network hosts one unit. The 
state of each unit is described by a dynamical variable 
b, called barrier. In each step, the unit with minimum 
b is mutated by updating the barrier. The effect of the 
mutation on the environment is taken into account by 
changing b also at all neighbors of the minimum site. 

General feature of ED models is the avalanche dynam- 
ics. The forward A-avalanches are defined as follows [Q. 
For fixed A we define active sites as those having bar- 
rier 6 < A. Appearance of one active site can lead to 
avalanche-like proliferation of active sites in successive 
time steps. The avalanche stops, when all active sites 
disappear again. There is a value of A, for which the 
probability distribution of avalanche sizes obeys a power 
law without any parameter tuning, so that the ED mod- 
els are classified as a subgroup of self-organized critical 
models 

The BS model was originally devised in order to ex- 
plain the intermittent structure of the extinction events 
seen in the fossil record In various versions of the 

BS model it was found, that the avalanche exponent is 
1 < r < 3/2, where the maximum value 3/2 holds in 
the mean-field universality class. On the other hand, in 
experimental data for the distribution of extinction sizes 
higher values of the exponent, typically around r 2± 2 are 
found ||T^. The avalanche exponent close to 2 was also 
measured in ricepile experiments [ p^ . While there are 
several models of different kind, which give generic value 



T = 2 ||l^.0, we are not aware of any ED model with 
such a big value of the exponent. 

The universality class a particular model belongs to, 
depends on the topology of the network on which the 
units are located. Usually, regular hypercubic networks 
§ or Cayley trees g are investigated. For random 
neighbor networks, mean-field solution was found to be 
exact Also the tree models |11 were found to be- 

long to the mean-field universality class. Recently, BS 
model on random networks, produced by bond percola- 
tion on fully connected lattice, was studied |l|]. Two 
universality classes were found. Above the percolation 
threshold, the system belongs to the mean-field univer- 
sality class, while exactly at the percolation threshold, 
the avalanche exponent is different. A dynamics chang- 
ing the topology in order to drive the network to critical 
connectivity was suggested. Similar model was investi- 



gated recently |20| in the context of autocatalytic chem- 
ical reactions. 

We present here a further step towards reality. In fact, 
one can find real systems in which not only the topol- 
ogy of connections evolves, but also the number of units 
changes. The network develops new connections, when a 
new unit is inserted, and if a unit is removed, its links are 
broken. This is a typical situation in natural ecologies, 
where each extinction and speciation event changes also 
the topology of the ecological network. The same may ap- 
ply to economics and other areas, where the range of in- 
teraction is not determined by physical Euclidean space. 
This problem was already partially investigated within 
mean-field BS model and also in several other models 
devised for description of biological evolution ]r7| , p2| -p4[ , 
which use different approaches than the extremal dynam- 
ics. 

The purpose of this Letter is twofold. First, to study 
the evolution of topology in a ED model with variable 
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number of units. Second, to demonstrate that the large 
value of the avalanche exponent can be observed if the 
topology of the underlying network evolves dynamically. 

We consider a system composed of varying number 
of units connected in a network. In the context of bio- 
logical evolution, these units are species. The dynamical 
rules of our model are the following. 

(i) Each unit has a barrier b against mutations. The 
unit with minimum b is mutated. 

(ii) The barrier of the mutated unit is replaced by a 
new random value b' . Also the barriers of all its neighbors 
are replaced by new random numbers. If b' is larger than 
barriers of all its neighbors, the unit gives birth to a 
new unit (speciation). If b' is lower than barriers of all 
neighbors, the unit dies out (extinction). As a boundary 
condition, we use the following exception: if the network 
consists of a single isolated unit only, it never dies out. 

This rule is motivated by the following considera- 
tions. We assume, that well-adapted units proliferate 
more rapidly and chance for speciation is bigger. How- 
ever, if the local biodiversity, measured by connectivity 
of the unit, is bigger, there are fewer empty ecological 
niches and the probability of speciation is lower. On the 
other hand, poorly adapted units are more vulnerable 
to extinction, but at the same time larger biodiversity 
(larger connectivity) may favor the survival. Our rule 
corresponds well to these assumptions: speciation occurs 
preferably at units with high barrier and surrounded by 
fewer neighbors, extinction is more frequent at units with 
lower barrier and lower connectivity. Moreover, we sup- 
pose that a unit completely isolated from the rest of the 
ecosystem has very low chance to survive. This leads to 
the following rule. 

(in) If a unit dies out, all its neighbors which are not 
connected to any other unit also die out. We call this 
kind of extinctions singular extinctions. 

From the rule (ii) alone follows equal probability of 
adding and removing a unit, while the rule (Hi) enhances 
the probabilty of the removal. As a result, the proba- 
bility of speciation is slightly lower than the probability 
of extinction. The degree of disequilibrium between the 
two depends on the topology of the network at the mo- 
ment and can be quantified by the frequency of singular 
extinctions. The number of units perform a biased 
random walk with reflecting boundary at riu = 1. The 
bias towards small values is not constant, though, but 
fluctuates as well. 

(iv) Extinction means, that the unit is removed with- 
out any substitution and all links it has, are broken. Spe- 
ciation means, that a new unit is added into the system, 
with a random barrier. The new unit is connected to all 
neighbors of the mutated unit: all links of the "mother" 
unit are inherited by the "daughter" unit. This rule re- 
flects the fact that the new unit is to a certain extent 
a copy of the original, so the relations to the environ- 
ment will be initially similar to the ones the old unit has. 



Moreover, if a unit which speciates has only one neigh- 
bor, a link between "mother" and "daughter" is also es- 
tablished. 
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FIG. 1. An example of network configurations in several 
successive time steps, from step 369 to 377. Units are repre- 
sented by small filled circles. The lines represent links between 
the units. The unit with the minimal barrier b is indicated 
by a filled square. 

The above described rules are illustrated in Fig. |^. 
Speciation occurs in the transition from step 369 to 370, 
371 to 372, 373 to 374 and 375 to 376. In the transitions 
from 372 to 373, 374 to 375 and 376 to 377 extinction 
occurs (last two include also singular extinctions). We 
can see that after speciation the neighbors of the new 
unit have one neighbor more than before, so that if 
increases, also the connectivity of the network grows. 

We investigated the evolution of the network by mea- 
suring time dependence of several quantities. We start 
the simulation with the initial condition — I. A typ- 
ical result is shown in Fig. where we show the time 
dependence of the number of units n^, average connec- 
tivity c, and the frequency of singular extinctions /g. The 
network can be split into disconnected clusters, as is illus- 
trated in the Fig. |l|. In Fig. || we show also the evolution 
of number of clusters in three size categories. We denote 
by ni the number of the smallest clusters, of size 2, by 
712 the number of medium-size clusters, larger than 2 and 
smaller or equal to nu/2, and by the number of clus- 
ters larger than The value — 1 means that most 
of the system is concentrated in a single cluster. 

We can see that singular extinctions occur in bursts. 
The passages without singular extinctions, where number 
of units evolves like a random walk (due to equal proba- 
bility of increase and decrease of number of units by 1) 
are interrupted by short periods, where rzu falls to small 
values and singular extinctions are intense. We can see 
that very often three events coincide: high frequency of 
singular extinctions, large number of clusters, especially 
of size 2, and the fact, that the largest cluster does not 
contain most of the network. We observed, that the mu- 
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tation occurs nearly all the time in the largest cluster. A 
similar effect was reported also in the Cayley tree models 
|pj| : the small isolated portions of the network are very 
stable and nearly untouched by the evolution. 
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FIG. 2. Time evolution of quantities describing the size and 
topology of the network, in a typical run. From top to bottom: 
the number of units riu (upper line) and average connectivity 
c (lower line), number of singular extinctions during 10^ steps, 
number of clusters of size 2, number of clusters of size larger 
than 2 and < number of clusters of size larger than 

nu/2. 

Figure || suggests that the number of units exhibits in- 
termittent drops. This corresponds qualitatively to punc- 
tuated equilibria seen in the fossil data. This feature 
is new in our model, when compared with previous ap- 
proaches within the BS model ||2^ but resembles the in- 
termittent features of models based on neural networks 
[17|, Lotka-Volterra equations and coherent noise 

[22|. In order to check this property quantitatively, we 
plot in Fig. ^ the distribution of changes in number of 
units during time interval At = 3 • 10^ steps. We can see 
that the distribution of drops (Ariu < 0) has a power-law 
tail, which confirms the intermittency. 

The distribution of forward A-avalanches is shown in 
Fig. ^. We found power-law distribution for Ac — 0.016 
with the exponent r = 1.98 ± 0.04. The value of the 
exponent was found by fitting the data in the interval 
(500,5 • 10^) H. Contrary to the BS and related mod- 
els, we found power law distribution with a different ex- 
ponent r' = 1.65 ± 0.05 also for A larger than about 0.4. 
(More precisely, we obtained the value of the exponent 



by fitting the data in the interval (500, 10*^) for A = 0.4 
and A = 0.6. Both values of A give the same result). The 
data suggest that for A > Ac the avalanche size distribu- 
tion exhibits two regimes, with crossover around certain 
avalanche size Scross- For small avalanches, s < Scross, 
the distribution is power law with exponent r, while for 
larger avalanches, s > Scross, power law with exponent r' 
holds. The crossover Scross grows when A approaches to 
Ac. 
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FIG. 3. Distribution of drops in number of units during 
the interval At = 3 ■ 10* steps. The full line is a power with 
exponent -3. 

The existence of avalanches for A close to 1 is related 
to the fluctuation of the number of units. We observed, 
that such avalanches start and end mostly when number 
of units is very small. Between these events the evolu- 
tion of the number of units is essentially a random walk, 
because singular extinctions are rare. This fact can ex- 
plain, why the exponent t' is not too far from the value 
3/2 corresponding to the distribution of first returns to 
the origin for the random walk. The difference is proba- 
bly due to the presence of singular extinctions. 




FIG. 4. Distribution of forward avalanche sizes, averaged 
over 12 independent runs. Each run lasted 3 • 10* steps. The 
values of A are (from bottom to top) 0.013, 0.016, 0.02, 0.03, 
0.05, 0.1, 0.2, 0.4, 0.6. The superscript > indicates that we 
count all events larger than s. 

The fact, that the number of units change, enables us 
to define the extinction size in a more realistic way than 
in the previous variants of the BS model. For fixed A we 
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count number of units which were present at the begin- 
ning of the A-avalanche and are no more present when 
the avalanche stops. This quantity corresponds better to 
the term "extinction size" used by paleontologists than 
the number of units affected by mutations, as it is defined 
in the BS model. 
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FIG. 5. Distribution of extinction sizes for A — 0.016, aver- 
aged over 12 independent runs. Each run lasted 3 • 10* steps. 
The full line corresponds to exponent Text = 2.32. The super- 
script > indicates that we count all events larger than s. 

For X ~ Xc — 0.016 the distribution of extinction sizes 
follows the power law with exponent Text = 2.32 ± 0.05, 
as is shown in Fig. |[ (The value of the exponent was ob- 
tained by fitting the data in the interval (10, 1000).) This 
value is larger than the exponent 2 observed in the statis- 
tics of real biological extinctions, but still closer than the 
values found in previous modifications of the BS model. 

The fact, that the network evolves and the number of 
units fluctuates leads to significantly larger values of the 
exponent than in the BS model, even greater that the 
experimental one, while variants of the BS model have 
values lower than the experimental one. This suggests, 
that the freedom in changing the topology in our model 
is exaggerated and in order to obtain a more realistic 
model of the biological evolution, we should look for some 
principles, which imply freezing of the topology, while 
allowing the species to be replaced by new ones. 

We studied several other modifications of our model, in 
order to check its robustness. For example, the link be- 
tween "mother" and "daughter" unit was established, or 
only certain fraction of the links connecting "mother" to 
its neighbors was inherited by "daughter" . These modi- 
fications affected some aspects of the network dynamics, 
but the avalanche and extinction statistics was not sig- 
nificantly different. 

To sum up, we formulated and studied the extremal 
dynamics model derived from the Bak-Sneppen model, 
which exhibits forward-avalanche exponent close to two, 
due to the annealed topology of the network. The extinc- 
tion size was defined in more realistic manner compared 
to previous approaches within the BS model and the ex- 
tinction statistics was found to obey a power law with 
exponent somewhat larger than two. The value found is 



closer to paleontological data than in the previous vari- 
ants of the BS model. 
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